A possible mechanism for cold denaturation of proteins 

at high pressure 
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We study cold denaturation of proteins at high pressures. Using multicanonical Monte Carlo 
simulations of a model protein in a water bath, we investigate the effect of water density fluctuations 
on protein stability. We find that above the pressure where water freezes to the dense ice phase 
(~ 2 kbar), the mechanism for cold denaturation with decreasing temperature is the loss of local 
low-density water structure. We find our results in agreement with data of bovine pancreatic 
ribonuclease A. 



Some proteins become thermodynamically unstable at 
low temperatures, a phenomenon called cold denatura- 
tion [H El This phenomenon has been mainly ob- 
served at high pressures, ranging from approximately 200 
MPa up to 700 MPa 0. An explanation of the P - T 
phase diagram of a protein with cold denaturation has 
been proposed [5j, but a microscopic understanding of 
the mechanisms leading to cold denaturation has yet to 
be developed, due in part to the complexity of protein- 
solvent interactions. 

Existing theories of folding and unfolding of diluted 
proteins consider hydrgphob icity as the driving force of 
protein stability S El- in the case of apolar 
macromolecules, hydrophobicity has been identified with 
the assembly and segregation of the macromolecule to 
minimize the disruption of hydrogen bonds among water 
molecules [(J [ijj, lll| . Water tends to be removed from the 
surface of apolar molecules, forming a cage composed of 
highly organized water molecules around the molecule, 
where the disruption of hydrogen bonds is minimized 
[l2|. The simplest hydrophobic model features an effec- 
tive attraction between hydrophobic molecules |l3|| . but 
does not reproduce cold denaturation. In order to obtain 
cold denaturation with this model, new studies 0, 0] 
had to insert a temperature-dependent attraction derived 
from experimental observations at ambient pressure . 
An explicit account of water around the hydrophobic 
molecules has also been considered in order to under- 
stand the cold denaturation process with temperature- 
independent interactions. Theoretical attempts modeled 
the effective water-protein interactions with the free en- 
ergy cost of excluding the solvent around the nonpolar 
molecule 0, . Numerical simulations based on these 
attempts have been applied to study the pressure denat- 
uration found in proteins 0. 

Not until recently has cold denaturation been studied 
at the microscopic level. Simplified models , based on 
a bimodal description of the ener gy o f water in the shell 
around the hydrophobic molecule , predicted cold de- 
naturation. Similar results were obtained using a lattice 
model of a random hydrophobic-hydrophilic heteropoly- 



mer interacting with the solvent [20j. Several models 
mimicking the interaction between water molecules and 
non-polar monomers have also been applied to the study 
of cold denaturation. [2l| . One possible reason for the in- 
ability of the previous models to capture both the molec- 
ular details of cold denaturation and the effect of pressure 
is the neglect of (i) correlations among water molecules 
near the freezing point, and (ii) the density anomaly 
due to the tetrahedral structure of the hydrogen bonded 
network. Here, we implement a two-dimensional lattice 
model of water that captures the above mentioned water 
properties • In the model, the possible orientations of 
water molecules are set by the allowed values of a g-state 
Potts variable Cj. Only when two neighbor molecules 
are in the correct orientation (<x; = aj) does a hy- 
drogen bond (HB), that increases the volume of the sys- 
tem by AV, form. This interaction mimics the increment 
of volume due to the incipient formation of a tetrahedral 
structure. If the two neighbor molecules are not 

in the correct orientation, the interaction of the particles 
does not imply any increment in volume. The Hamil- 
tonian for our model of water-water interaction may be 
written as 221 
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were J > is the scale of the interaction between water 
molecules upon tetrahedral network formation. The to- 
tal volume of the system is given by V = Vq + NhbAV, 



where Nhb = 



(id) 



7j is the total number of hydro- 



gen bonds with AV > in the system. The sum 
extends only to nearest neighbors, implying that two wa- 
ter molecules cannot form a hydrogen bond with Ay > 
if they are separated by one residue of the protein. The 
enthalpy of the system (Hhb + PV) is given by 

H HB + PV = - ( J - PAV) J2 W i > (2) 

(id) 

where P is the pressure applied to the system. 
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Our model solvent features a limiting pressure P c = 
J/AV. Above P c , we find that Nhb decreases as we 
decrease the temperature, and the water model under- 
goes a transition to a state where all hydrogen bonds 
with AV > are broken. Below P c , we find that Nhb 
increases as we decrease the temperature, and the wa- 
ter model undergoes a sharp transition, at T — T c — 
(J - PAF)/(ln(l + y^)) jig, to a state where all hy- 
drogen bonds with AV > are formed. Thus, our water 
model reproduces the freezing of water to low and high 
density ice, since for P < P c (P c « 200MPa in real water) 
water freezes to the low density ice Ih, and for P > P c , 
water freezes to the high density ice II [23|. A relation 
between these two phases of ice and protein folding has 
already been suggested from a thermodynamic point of 
view |24j . 

We model the protein as a self-avoiding random walk 
embedded in a water bath. For simplicity, we consider 
a non-polar homopolymer that interacts with water via 
the partial ordering of water molecules, forming hydrogen 
bonded structures around the protein. We mimic the 
interaction using the Hamiltonian, 



Hp = J r nHB(n r . 



^2 niU ^ ' 
{id) 
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where the parameter J r > is the strength of the repul- 
sive interaction and uhb = Nhb /N wate r is the number 
density of hydrogen bonds with AV > 0, N wate r is the 
number of water molecules. The water-protein repulsion 
increases as the water molecules tend to form the tetra- 
hedral, low-density, hydrogen bonded network, where an 
unfolded apolar macromolecule is unlikely to be embed- 
ded, i^max is the maximum number of residue— residue 
contacts and J^V ^ riiUj is the number of residue-residue 
contacts, where rii = 1 if the lattice position i is oc- 
cupied with a residue, and zero otherwise. Therefore, 
n m ax — n i n j is a measure of protein compactness, 

and equals zero when the protein is maximally compact. 
Thus, Eq. states that the hydrophobic repulsion driv- 
ing the protein to a compact state is equal to zero when 
the protein is maximally compact Q^<ij) n i n j ~ n max ) 
or when the water forms the high-density bond network 
(n H B ~ 0). 

We hypothesize that the inability of water molecules 
to arrange in the low density ice-like structures is the 
principal mechanism responsible for protein cold denatu- 
ration. At low pressures (P < P c ) and low temperatures, 
water molecules form a low density hydrogen bonded net- 
work, so the protein is forced to adopt a compact state. 
At high pressures (P > P c ), the water is not able to form 
the low density network and forms a more dense state. 
In this case the effective repulsion between the residues 
and the solvent decreases, and water molecules penetrate 
into the protein core, unfolding the compact state. Our 
hypothesis is supported by the experimental observations 



[J] that cold denaturation exists mainly at high pressures 
(of the order of kbars), where water only freezes in the 
dense ice II phase [25^ . 

Next we demonstrate that our model of a protein em- 
bedded in a water network with an enthalpy given by 



W = H p +Hhb +PV 



(4) 



gives rise to both cold and warm denaturation of the 
protein and agrees with experimentally-observed protein 
denaturation at high pressures. Since the energy land- 
scape of the protein interacting with the water network 
is characterized by a multitude of local minima, we per- 
form multicanonical Monte Carlo simulations to avoid 
transient trapping of our water-protein system in local 
energy minima at low temperatures. Specifically, we use 
the multiple-range random walk algorithm |28j to cal- 
culate the density of states. We adopt the algorithm 
in order to embed the self-avoiding protein into the lat- 
tice, and to calculate the two-parameter density of states 
9(Nhb, N c ), where N c = i s the number of 

residue-residue contacts. From the density of states we 
calculate the temperature and pressure dependence of the 
average number of residue-residue contacts 



_ e -W(N HB .N c )/T 

Nc=J2 J2 N -9( N hb,N c ) , 

Nhb N 



(5) 



where Z is the partition function. We perform Monte 
Carlo simulations of a system of 383 water molecules 
and a protein consisting of 17 non-polar residues with 
periodic boundary conditions. Fig. ^ shows the depen- 
dence of N c /n max on temperature for different values 
of the pressure both above and below P c . The calcu- 
lated density of states g(NnB, N c ) converges to the true 
value with an accuracy of the order of 10 -5 . The value 
of N c ln m ax ranges from one (maximally compact pro- 
tein) to approximately 0.71 (which is the average number 
of residue-residue contacts found at high temperatures). 
Only when P > P c do we observe the cold denaturation 
of the protein. 

We also reconstruct the phase diagram of the water- 
protein system in the P — T plane (Fig. [2I> - We consider 
the protein to be in the collapsed state if 96% of all pos- 
sible contacts are formed, i.e., if N c /n max > 0.96. For 
each pressure value, the freezing lines of water shown in 
Fig. are given by the temperature at which we observe 
a maximum in the specific heat of the water bath. We 
compare our findings to experimental observations |29T | 
for bovine pancreatic ribonuclease A studied by 1 H NMR 
spectroscopy. We find remarkable qualitative agreement 
between the experimental and numerical P~T phase dia- 
grams. In both experimental and numerical P — T phase 
diagrams, we observe that cold denaturation occurs at 
high pressures and, as we lower the temperature, close 
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FIG. 1: Normalized number of residue-residue con- 
tacts vs. temperature for pressure values PVo/J = 
0, 0.25, 0.5, 0.75, 1, 1.1, 1.125, 1.15, 1.175, 1.2, 1.25, 1.4, 1.5, 2. 
A magnified region where the cold denaturation takes place 
is shown in the inset. The dotted line corresponds to 
Nc/rimax > 0.96. We represent the curve corresponding to 
P = P c = J I AV by a bold line. The values of the parameters 
used are: J = 1, J r = 10, AV = 1 and q = 10. 



FIG. 2: P-T phase diagram for the protein derived from Fig.0 
Dashed lines indicate the freezing lines for model water. Wa- 
ter freezes in low density ice Ih for PVo/J < 1 and in dense 
ice II for PVo/J > 1. In the inset we present the experi- 
mental results obtained by Zhang et al. |29| for the bovine 
pancreatic ribonuclease A. Two typical configurations of the 
protein are shown, one in the compact state and the other in 
the denaturated state. 



to the water-ice II freezing line and in the region where 
water molecules are not capable of forming low density 
ice-like structures. In addition to the study of ribonu- 
clease A [i^ , cold denaturation at very high pressures in 
the kbar range has also been observed in chymotrypsino- 
gen 0, myoglobin 0, staphylococcal nuclease [Ho| and 
has been proposed as the principal mechanism for the 
observed pressure-inactivation of bacteria, such as Es- 
cherichia coli |3l| . 

Not all proteins behave equally as we decrease tem- 
perature at high pressure. In particular, there are some 
proteins that do not exhibit cold denaturation. We re- 
produce the variability of protein dynamics at high pres- 
sure and low temperature by varying the hydrophobic 
parameter J r to lower values, effectively impeding a sta- 
ble compact state for pressures above the P = P c line. 
In Fig. |3| we present the phase diagrams obtained for dif- 
ferent values of J r , ranging from two to 20. The shape 
of the phase diagram changes as we increase the value of 
the repulsive interaction J r , allowing stabilization of the 
compact state and cold denaturation above the P = P c 
line. 

Within the framework of our model, we reproduce the 
experimentally-observed thermodynamics of cold denat- 
uration, but we cannot address the kinetics of this pro- 
cess. It could be also feasible to investigate in future work 
the dynamics with a more sophisticated model, where we 
consider not just the average number of hydrogen bonds 
with AV > 0, but the actual numbers for each water 



ICE II . 




FIG. 3: P-T Phase diagram for proteins with J r = 2, 5, 10, 20. 
Dashed lines indicate the computed freezing lines for water. 
Doted line indicates the P = P c = J/AV line. 

molecule that is a neighbor to a residue. 

Recent computer simulations studies | 32j with all-atom 
models have studied the effect of pressure and average 
density on the hydrophobic effect. The authors per- 
formed simulations of two Lennard- Jones (LJ) particles 
in the TIP4P water model At constant tempera- 

ture, they found that the aggregation of the two particles 
is favored with a moderate increase of pressure, or analo- 
gously, with a moderate increase of water density. How- 
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ever, this effect is reversed for pressures in the kbar range, 
so that aggregation becomes unstable with increasing wa- 
ter density. These studies support our results of a critical 
pressure above which an increase of pressure leads both 
to an increase of water density (because of the reduction 
of the number of hydrogen bonds with AV > 0) and 
to destabilization of the model protein. Following these 
results and based on previous studies [lfj, Shimizu et 
al. have shown that three-body interactions have a 
destabilization effect on the aggregation of three L J par- 
ticles. However, the inclusion of these interactions into 
a model of a more complex protein did not lead to sig- 
nificant changes. Finally, all-atom simulations recently 
addressed the pressure denaturation of proteins Hi,but 
they have not provided conclusive evidence. 

We conclude that the effect of pressure on water den- 
sity is key for understanding cold denaturation of pro- 
teins. The density anomaly of water arises from the low 
density hydrogen bonded structures responsible for the 
hydro pho bic effect, driving the protein to a compact state 
[a IsL Unl flijl . At extreme pressures above P c , lowering 
the temperature implies an increasing free energy cost 
to form a hydrogen bond with AV > 0, so the density 
anomaly disappears. In this scenario, the hydrophobic ef- 
fect decreases and cold denaturation occurs. Our model 
supports this mechanism. Also, a specific arrangement 
of amino acids in the protein structure, determined by 
amino acid interactions, dictates the dynamics of pro- 
teins at low temperature and high pressure, thus making 
some proteins more stable than others at these P — T 
conditions. 
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